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We have studied the efltect of resonant electronic state coupling on the formation of ultracold 
ground-state *^Rb2. Ultracold Rb2 molecules are formed by photoassociation (PA) to a coupled pair 
of 0^ states, 0j(Pi/2) and 0^(P3/2), in the region below the 5S + 5Pi/2 limit. Subsequent radiative 
decay produces high vibrational levels of the ground state, X ^E^. The population distribution of 
these X state vibrational levels is monitored by resonance-enhanced two-photon ionization through 
the 2 ^Tit, state. We find that the populations of vibrational levels ii"=112— 116 are far larger than 
can be accounted for by the Franck-Condon factors for 0^(Pi/2) — > X transitions with the 
0j(Pi/2) state treated as a single channel. Further, the ground-state molecule population exhibits 
oscillatory behavior as the PA laser is tuned through a succession of 0^ state vibrational levels. Both 
of these effects are explained by a new calculation of transition amplitudes that includes the resonant 
character of the spin-orbit coupling of the two Oj states. The resulting enhancement of more deeply 
bound ground-state molecule formation will be useful for future experiments on ultracold molecules. 

PACS numbers: 32.80.Pj, 32.80.Qk, 34.50.Rk. 



I. INTRODUCTION 



40 



30 



o 

20 



10 











~600nm 








55+4D 






55+5P 










SE 


PA 

55+55 





10 
R(A) 



15 



20 



FIG. 1: (color online) Selected potential curves for Rb2 from 
Ref. y^], showing our scheme for forming ultracold *^Rb2 
molecules in the X^E+ state. Photoassociation (PA) is used 
to excite high-u levels of the coupled states, 0^ (P1/2) and 
(P3/2), converging to the 5S'-|-5Pi/2 and 5S'-|-5P3/2 atomic 
limits. These states correlate at short range to the A ^Ej and 
h "^IIu states. Spontaneous emission (SE) from these levels 
populates the ground X state. Pulsed laser light at 600 nm 
is used for ionization detection. 



'Present address: Akamai Physics, Inc., Las Cruces, NM 88005 



In recent years, the technique of photoassociation (PA) 
has become a powerful method for the formation and 
study of ultracold molecules produced from ultracold 
atoms d i]. Although PA can form molecules only 
in electronically excited states, subsequent radiative de- 
cay can produce ground-state molecules, mainly in their 
highest vibrational levels. To enhance the formation of 
ground-state molecules in levels well below the dissoci- 
ation limit, it would be ideal to photoassociate to an 
excited state that has a bimodal vibrational probability 
distribution, with one peak at long range to enhance the 
photoassociation rate, and another at shorter range to 
enhance radiative decay into deeply-bound levels of the 
ground electronic state. Theoretical and experimental 
work on ultracold Cs2 [1, 0] has shown that such states 
can exist because of resonant spin-orbit coupling between 
the vibrational series of the two 0+ states that converge 
to the (65' -I- 6P1/2) and (GS' -I- 6P3/2) fine-structure hm- 
its. In Ref. 3 the formation of more deeply-bound 
vibrational levels of Cs2 is greatly enhanced when this 
coupling is at its maximum, and it was proposed that 
this approach might be equally applicable for other al- 
kali dimers. There have been several additional studies 
of the effects of the resonant coupling phenomenon on 
photoassociation spectroscopy and efficient formation of 
ultracold molecules d, H, 0, HI ■ 

In Rb2 the spin-orbit perturbations between the pair 
of 0+ states converging to (SS'-I- 5P1/2) and (55* -I- 5P3/2) 
have been studied in detail 0, [1] , and agreement between 
theory and experiment for the energy levels is very 
good. However, the effects on ground-state Rb2 forma- 
tion have not been investigated until quite recently. In 
2006, anomalous vibrational distributions were reported 



2 



that could not be fully explained using Franck-Condon 
factors calculated from adiabatic potentials A re- 
cently completed comparative study of molecule forma- 
tion in ^^Rb and ^^Rb attributes most of the differences 
between the isotopes to resonant coupling [10]. In this 
paper we report a detailed experimental and theoreti- 
cal study of resonant enhancement of ground-state ^^Rb2 
molecule formation, with an emphasis on the energy de- 
pendence for a range of vibrational levels of the 0+ state. 
As shown in Fig. [U we photoassociate pairs of colliding 
^^Rb atoms into high molecular vibrational levels of the 
two 0+ states in the region just below the 55' -I- 5Pi/2 
asymptote. Some of these molecules decay radiatively 
to the X^S+ state, forming either bound molecules in 
high-w" levels or unbound pairs of atoms. The bound 
ground-state molecules are detected state-selectively by 
two-photon ionization via the 2 state. 

We focus particularly on the variation of molecule for- 
mation rates for the vibrational levels v" = 112—116 
as a function of the vibrational level v' of the coupled 
0^ states to which PA takes place. A coupled-channels 
model is used to predict Franck-Condon factors and tran- 
sition moments for the 0^ ^ X ^ radiative decay pro- 
cess, as well as the detection process. These results are 
compared with measurements of the relative formation 
rates after carefully normalizing them to remove depen- 
dencies on the PA rates and extraneous experimental fac- 
tors. 



II. THEORY 

The coupled 0+(Pi/2) and 0+(P3/2) states arise due to 
spin-orbit coupling of the short-range j4^E+ and b^Iiu 
states. Dipole-allowed radiative decay to the ground 
Ar^I]+ state is possible only for the portion of the 0+ 
wavefunction having character. Vibrational eigen- 

functions of the 0+ excited states, the X ^S]+ electronic 
ground state and the 2 excited state are obtained in 
order to calculate Franck-Condon factors needed to pre- 
dict the relative efficiency of formation and detection of 
ground-state molecules in various vibrational levels v" . 



A. Model for the molecule formation process 

In the diabatic basis, the Hamiltonian describing nu- 
clear motion in the electronically excited Oj(Pi/2) and 
0« (-P3/2) states reads 

V <o^(^) t + H3n„(i?)-W^|^(i?)y ' 

where R denotes the internuclear coordinate, T the corre- 
sponding nuclear kinetic energy operator, Vi{K) the re- 
spective potentials, and WIq{R) the diagonal and off- 



diagonal spin-orbit coupling functions. The hyperfine in- 
teraction in the excited state can be neglected because 
PA detunings of several wavenumbers are considered. We 
are interested in the formation of ground-state molecules 
with binding energies E^jh > 3cm^^, much larger than 
the ground-state hyperfine splitting of 0.1 cm^^ for *^Rb. 
Therefore, the hyperfine interaction in the ground state 
can also be neglected, and the A^I]+ ground state is 
hence described in a single channel picture. The 2^Ej 
state can also be treated as a single channel. The corre- 
sponding Hamiltonians are 

H, = t + y,(p) (2) 

with i = A:1I]+,21I]+. 

The potentials Vi{R) for the and 6^n„ states 

and the spin-orbit coupling functions WI^q{R) are taken 
from Ref. [7| . The ground-state potential is obtained by 
connecting the short range potential curve of Ref. ll] to 
the long-range expression Cq/R^ + Cs/R^ + Cio/R^° with 
the Cj coefficients found in Ref. ^1^]. The repulsive wall 
of the ground-state potential is adjusted to yield a singlet 
scattering length of 2400 ao for *^^Rb [l3|. Thepotential 
curve for the 2 ^Ej state is taken from Ref. [ij]. 

The Hamiltonians, Eqs. (fT]|2l) . are represented on a co- 
ordinate space grid employing a mapped grid method [isl . 
Il6i] . Vibrational energy levels and wavefunctions are ob- 
tained by diagonalization, and Franck-Condon factors 
are calculated by evaluating the appropriate integrals. 
For the 0+ AT^E^ transitions, Franck-Condon fac- 
tors were compared to transition dipole matrix elements 
where the transition dipole was included in the integrals 
as a function of internuclear coordinate • No signifi- 
cant differences were found in the relative transition am- 
plitudes, consistent with the expected nearly constant 
long-range behavior of the transition dipole for the range 
of PA detunings considered here. 



B. Resonant enhancement of ground-state molecule 
formation with binding energies of 2 — 20 cm~^ 

It is well-known [1, [11, ^ that the 0+(Pi/2) and 
0+(P3/2) states of Rb2 constitute a classic example of 
strongly coupled electronic states. This implies that 
throughout the vibrational spectrum, some of the 
eigenfunctions are strongly perturbed by the resonant 
coupling while others show fairly regular behavior. Fig. [2] 
compares strongly perturbed and regular vibrational 
eigenfunctions with binding energies of 9.39 cm^^ and 
11.27 cm~^, respectively. 

The spectral positions at which the two 0^ components 
are most strongly coupled can be readily identified, e.g. 
by inspection of the rotational constants B{v) (cf. Fig. 6 
of Ref. 0). The behavior of the square of the wave- 
function of Fig. [5] (a) arises from the resonant coupling 
of the eigenfunctions of the two 0„ potentials, showing 
two classical turning points for each of the two potentials. 
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FIG. 2: (color online) The squares of the vibrational eigen- 
functions of the 0^ coupled states for binding energies of 
9.39 cm~^ and 11.27 cm~^ (red solid line: triplet component, 
blue dotted line: singlet component). The vertical position 
of the square of each wavefunction corresponds to its binding 
energy. The potentials (thick solid lines) are drawn in order 
to illustrate how the structure of the coupled wavefunctions 
can be understood in terms of classical turning points. 



The enhanced probability ampHtude at the outer classical 
turning point of the upper 0^(^3/2) potential, at about 
R = 23.5 oq, is expected to lead to large Franck-Condon 
factors with vibrational levels in the electronic ground 
state. The variation of the square of the wavefunction 
in the vicinity of this classical turning point is empha- 
sized in the inset of Fig. [5] (a). By contrast, the square of 
the vibrational wavefunction in Fig. [2] (b) shows almost 
regular behavior with only a small perturbation visible 
around the outer classical turning point of the 0+(P3/2) 
potential. Note that the scale of the inset of Fig. [2] (a) is 
three times the scale of the inset of Fig. [5] (b). 

The Franck-Condon factors for each of the two 0+ 
eigenfunctions of Fig. [5] with the highest vibrational lev- 
els of the X^'E^ state are presented in Fig. [3l For reg- 
ular vibrational wavefunctions, one would expect larger 
0^ binding energies (corresponding to larger PA detun- 
ings) to result in larger Franck-Condon factors with the 
ground-state vibrational levels The Franck-Condon 
factor for the second to last X^T,'^ level {v" = 122) is in- 
deed slightly larger for the 0+ level bound by 11.27 cm^^ 
than for the level at 9.39 cm^^. The probability density 
of the level v" = 122 has its outermost peak at about 
i? = 60 Oq, just like the 0+ level with Sbind = 11.27 
cm"^ (cf. Fig. [2]), while the outermost peak of the 0+ 
level with i?bind = 9.39 cm~^ occurs at larger distances, 
about i? = 64 flQ. The main contribution to the Franck- 
Condon factors with v" = 122 can therefore be attributed 
to the long-range portions of the wavefunctions. These 
long-range maxima are also responsible for efficient pho- 
toassociation. 

In addition to this long-range contribution to the 
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FIG. 3: (color online) Franck-Condon factors between 
ground-state vibrational levels w" and the 0^^ vibrational lev- 
els with binding energies of 9.39 cm~^ (top) and 11.27 cm~^ 
(bottom). Enhancement of Franck-Condon factors by an or- 
der of magnitude due to the resonant coupling is observed for 
the X^T.+ levels v" = 113 (Sbind = 6.9 cm"^) and v" = 114 
(i^bind = 5.0 cm-i). 



Franck-Condon factors, the enhanced probability den- 
sity at intermediate distances (about R = 23.5 oq in 
Fig. ^ due to the resonant coupling leads to signifi- 
cant Franck-Condon factors with the ground-state levels 
v" = 112-119, in particular with v" = 113 and 114. 
Since the 0^ eigenfunction with Sbind — 9.39 cm^^ is 
strongly perturbed by the resonant coupling, its Franck- 
Condon factors with v" = 113 and 114 are larger by one 
order of magnitude than those of the 0^ eigenfunction 
with ii^bind — 11.27 cm^^. The enhancement of Franck- 
Condon factors for the levels v" = 112 — 119 due to strong 
perturbation of some of the 0„ vibrational wavefunctions 
will result in efficient formation of ground-state molecules 
in these levels which can be probed by varying the PA 
detuning. This is the first important conclusion to be 
drawn from Figs. [2] and [31 The effect is obvious in the 
solid curve and red diamonds of Fig. |31 which shows the 
integrated Franck-Condon factors from each 0+ vibra- 
tional level to ground-state levels with v" = 112 — 116 
(the levels v" = 117 — 119 are omitted because they can- 
not be efficiently detected in our experiment for reasons 
explained below) . Here the integration over a small range 
of v" smooths out the nodal oscillations that would be 
seen in the Franck-Condon factors to any single ground- 
state level v" , leaving only the overall structure due to 
resonant coupling of the 0,^ states. 

Moreover, as seen in Fig. [3l the Franck-Condon factors 
of the coupled-state 0^ eigenfunction of binding energy 
-Bbind = 9.39 cm~^ with the ground state levels v" = 113 
and 114 differ only by a factor of about 2.7 and 2.3, re- 
spectively, from the Franck-Condon factor with the sec- 
ond to last level, v" = 122. While the levels v" = 113 and 
114 are highly excited, their binding energies are much 
larger than that of v" = 122 (i?bind = 9.16 cm~^ and 
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FIG. 4: (color online) Solid curve with red diamonds: Sums 
of calculated Franck-Condon factors from each of 25 vibra- 
tional levels of the coupled states to ground-state levels 
with v" — 112 — 116, revealing strong oscillations in ground- 
state molecule formation. The Ot level energies are shown as 
a function of PA laser detuning below the 5S + 5Pi/2 limit for 
comparison with experimental data. Crosses: Experimental 
data set 1. Circles: Experimental data set 2, taken on a differ- 
ent day under improved conditions as explained in Section ITVl 
The typical systematic error bar shown is based on upper and 
lower bounds for systematic corrections as explained in the 
text. 

6.86 cm~^, respectively, compared to « 10~^cm~^ 
The outermost peaks of the corresponding vibrational 
wavefunctions occur at about 21.4 oq and 22.2 gq. This 
leads to the second important conclusion: Molecule for- 
mation in these more strongly bound ground-state levels 
is almost as efficient as in the very last levels just be- 
low the dissociation threshold. This appears to resolve a 
puzzle raised in Ref. [§], where it was observed that the 
vibrational populations of levels near v" = 112 are much 
higher than can be explained by Franck-Condon factors 
calculated using the uncoupled 0+(Pi/2) potential curve. 

III. EXPERIMENT 

The experimental setup has previously been described 
in detail in Ref. Q. It consists of a ^^Rb MOT, which 
is loaded to a density of approximately 10^^ cm~^ using 
a "dark-SPOT" ^ configuration. The estimated MOT 
temperature is 200 /iK. A cw tunable Ti:Sapphire laser 
(Coherent 899-29) with an output power of more than 400 
mW is focused into the MOT for photoassociation. The 
photoassociated ultracold molecules spontaneously decay 
to the ground state, X which is subsequently de- 

tected by resonance-enhanced two-photon ionization us- 
ing a pulsed dye laser (Continuum ND 6000) pumped by 
a 532 nm Nd:YAG laser at a 10 Hz repetition rate. The 
detection laser beam has a 1 mm diameter and a typical 
pulse energy of 3 mJ with a pulse duration of 7 ns. The 
molecular ions that it forms are accelerated into a chan- 
neltron detector and discriminated from atomic ions by 
their time of flight. 



With the PA laser frequency fixed on a resonance, the 
pulsed-detection laser can be scanned in the range of 601- 
603 nm to obtain the vibrational spectrum of the ground- 
state molecules. Since the theory of resonant coupling 
predicts an oscillatory behavior of the number of formed 
molecules as a function of PA laser detuning, we repeat 
the experiment for nine different PA laser detunings, each 
time obtaining a detection laser scan. To compare more 
clearly the effect of resonant coupling for these different 
PA laser detunings, we normalize each detection spec- 
trum by the fractional PA trap loss measured at the same 
PA laser detuning. More details regarding the normal- 
ization are given in Section IV. 

A limitation of this detection method is that it can- 
not detect the very highest vibrational levels, because 
they can be photodissociated indirectly via off-resonant 
re-excitation by the PA laser after they are formed. This 
indirect photodissociation process is discussed in detail 
in Ref. 9]. There it was observed that only levels up 
to v" = 117 can be observed, even though the formation 
rate for high-i;" levels should have been very large. In 
the present work we confine our analysis to levels with 
v" < 116, for which we estimate that losses due to this 
process should be minor. 

Another experimental limitation is that a small back- 
ground signal level is observed even in the absence of the 
PA laser, probably due to formation of high-u" ground- 
state molecules by the MOT laser itself. After noting 
this issue in our first data set, we included in a second 
data set explicit measurements of this background signal, 
taken by scanning the detection laser without the pres- 
ence of the PA laser. This background level can be sub- 
tracted from each ion spectrum. However, we note that 
the first data set cannot be so corrected, and also that 
the subtraction may over-correct for the real background 
contribution, since the presence of the PA laser may cause 
indirect photodissociation of the high-w" molecules com- 
prising the background, causing it to diminish whenever 
the PA laser is present. Thus we use in our analysis the 
average of the uncorrected and corrected signals, and we 
take the difference between the two to be upper and lower 
bounds when estimating the experimental uncertainties. 



IV. RESULTS AND ANALYSIS 

Two sets of experimental data were taken on different 
dates. Detection laser spectra were obtained for each of 
nine different detunings ApA of the PA laser, including 
the four examples shown in Fig. O The data include all 
of the PA transitions in the range ApA = 6—17 cm^^ that 
are sufficiently strong and free from overlaps with unre- 
lated transitions to allow reliable results. In this range, 
theory predicts two pairs of maxima and minima for for- 
mation of molecules due to the resonant coupling phe- 
nomenon. Assignment of the vibrational levels ?;"=112- 
116 is based on our previous study of the formation of 
the singlet molecules in the ground X^T,'^ state @. 
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FIG. 5: (color online) Comparison of experimental detection 
laser scans for PA detunings corresponding to the minimum 
and maximum effects of resonant coupling as predicted from 
the calculated curve in Fig. |3] Ionization signals for bind- 
ing energies of (a) 6.02 cm"^ (Min) (b) 9.39 cm"^ (Max) (c) 
11.27 cm"^ (Min), and (d)16.27cm"^ (Max) are shown. For 
easy comparison, the vibrational levels v" are marked by ver- 
tical lines. 



To show the enhancement of ultracold molecule forma- 
tion in the X ground state due to resonant coupUng, 
we compare the production of ground-state molecules for 
different PA laser detunings. However, the production 
rate involves both the efficiency of the PA step and the 
efficiency of the spontaneous emission (SE) step, as is 
apparent in Fig. [T] Only the latter process is of in- 
terest to the present analysis. Fortunately the effects 
of the two processes can be distinguished by using trap 
loss measurements to determine the PA rate. Produc- 
tion of excited molecules by PA results in trap loss with 
nearly 100% probability, either by formation of ground- 
state molecules or by decay back into the continuum 
with enough kinetic energy gain to escape from the trap. 
"When the PA laser is tuned to a resonance, the formation 
rate R of excited-state molecules is given by 



2 ' 



(3) 



where L is the MOT loading rate, is the steady- 
state number of atoms in the MOT in the absence of 
the PA laser, and AA^ is the change in number due to 
the PA laser. Therefore, the formation rate of excited- 
state molecules is proportional to the fractional trap loss 
with the assumption that the MOT parameters are fairly 
constant. Thus we can correct for variation in the PA 
rate at different PA laser detunings, by normalizing the 
ground-state molecule signal by the corresponding frac- 
tional trap loss. In a typical data set, this fractional loss 
varies from 18% to 27% over our range of detunings. As 



explained previously, we also sum together the signals 
(peak heights) from ground-state levels with w"=112-116 
to average over nodal fluctuations in the Franck-Condon 
factors. 

A small systematic error can occur because the detec- 
tion laser beam intensity is not always sufficient to sat- 
urate the first step of the detection process, the bound- 
bound X 2 transition. Therefore, we con- 
sider two limiting cases. First, we assume that the tran- 
sition is completely saturated. Then the ionization signal 
should directly reflect the Franck-Condon factors for the 
0+ ^^^g" transition. Second, we assume that there 
is no saturation at all. In this case the ionization signal 
must be normalized by the Franck-Condon factors for 
the flrst detection step, X — > 2 ^Tj'^. Because the 
extent of saturation was not well-controlled during the 
measurements, we use the average of these two limiting 
cases as our final experimental result. The difference be- 
tween them, which fortunately is small compared to the 
random experimental scatter, is taken as the associated 
systematic uncertainty. 

In Fig. [4l the final experimental results for radiative 
decay rates to levels with v" = 112 — 116 are compared 
with theory. Because the overall scale is arbitrary, we 
choose for each data set a multiplicative scaling factor se- 
lected to minimize the standard deviation between theory 
and experiment. The typical systematic error bars shown 
in Fig. 3] is estimated by taking the quadrature sum of 
the systematic uncertainties due to background subtrac- 
tion and incomplete detection laser saturation. It does 
not include the effects of random variations and day-to- 
day variations in the MOT, which are apparent in the 
scatter between the two independent data sets. 

The overall oscillatory behavior due to resonant cou- 
pling is very apparent in Fig. |4l and the positions of the 
minima and maxima coincide closely with the theoret- 
ical predictions. However, the ratio between the max- 
imum and minimum molecule formation rates is some- 
what smaller than predicted by theory, about 5:1 as op- 
posed to 20:1. This may arise from additional pertur- 
bations not taken into account in our two-state model, 
although it also may be due in part to experimental is- 
sues such as the effects of weak overlapping transitions 
and the saturation effects described above. "We also note 
that the relative populations of the individual vibrational 
levels v" exhibit variations different from the predicted 
Franck-Condon factors, even though the broader distri- 
butions are well-predicted, indicating that the model is 
not quite sufficient to predict the exact nodal positions 
in the vibrational wavefunctions. 



V. CONCLUSIONS 

We have presented clear theoretical and experimen- 
tal evidence for resonant enhancement in the photoas- 
sociative formation of ultracold ^^Rb2 molecules in the 
ground X^Yi^ state. The basic mechanism, enhance- 
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ment of radiative 0^ ^ X decay by resonant coupling of 
the 0+ (1/2) and 0+ (3/2) states, is similar to that pre- 
viously observed for Cs2 A range of PA transitions 
was investigated with detunings of 6 — 17 cm~^ below the 
55 + 5Pi/2 atomic asymptote. A mapped grid method 
was used to obtain vibrational wavefunctions to predict 
relative molecule formation rates, and experiments were 
performed to measure these rates. The oscillatory be- 
havior of the ground-state molecule population has been 
revealed more clearly than in the prior experimental work 
on Cs2. This was accomplished by normalizing the ex- 
perimental results to remove variations in PA rates and 
by integrating the formation rates for several vibrational 
levels to remove local nodal oscillations in the Franck- 
Condon factors. Agreement between theory and exper- 
iment is quite good, except that the observed contrast 
of the resonant enhancement phenomenon is somewhat 
less than predicted. The enhanced production of ground- 
state molecules in the vibrational levels v" = 112 — 116 



should be quite useful for future experimental studies. In 
particular, vibrational levels with binding energies of sev- 
eral wavenumbers present a much better starting point 
for additional Raman transitions to obtain molecules in 
their vibronic ground state than levels with sub-cm~^ 
binding energies [2l| . 
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